2, 2008 7:7 WSPC/INSTRUCTION FILE Lesmak'Meson2006v2 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



FINAL STATE INTERACTIONS IN B ^ nnK and B KKK 

DECAYS 



L. LESNIAK\ B. EL-BENNICH^, A. FURMAN^, R. KAMINSKI\ B. LOISEAU^ 

^ Division of Theoretical Physics, The Henryk Niewodniczanski Institute of Nuclear 
Physics, Polish Academy of Sciences, 31-342 Krakow, Poland 
Leonard. LesniakQifj. edu.pl 
Lahoratoire de Physique Nucleaire et de Hautes Energies (IN2P3-CNRS-Universites 
Paris 6 et 7), Groupe Theorie, 
Universite Pierre et Marie Curie, 4 PI- Jussieu, F-75252 Paris, France 
ul. Bronowicka 85/26, 30-091 Krakow, Poland 

Analysis of charged and neutral B meson decays into 'k'^tt~ K, K^K~ K and K^K^K^ 
is performed using a unitary representation of the tttt and KK final state interactions. 
Comparison of the theoretical model with the experimental data of the Belle and BaBar 
Collaborations indicates that charming penguin contributions are necessary to describe 
the B -> /o(980)/s: and B p(J7Q)° K decays. 
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1. Introduction 

We report on some studies of three-body B meson decays into tt-kK and KKK 
final states. In these reactions one can find an evidence of the direct CP violation 
similar to that recently discovered in two-body decays into t:^K^ . Loop- type 
weak decay diagrams known as penguin terms play an important role in these 
decays. In particular we study contributions of the charming penguin amplitudes 
responsible for long-distance effects present in the final state interactions. Strong 
interactions between pairs of pions and kaons in the final scalar-isoscalar state have 
been described in Ref. Here we extend this approach to the tttt interactions in 
P-wave. This enables us to obtain a unitarized description of the final state tttt 
interactions from threshold up to about 1.2 GeV. The opening of the KK threshold 
near 1 GeV is included in a natural way since both tttt and KK S-wave channels 
are coupled in our model. 

2. Decay amplitudes 

We take into account two components in the weak transition amplitudes governed 
by 6-quark decays into suu, sdd and sss, where u, d and s denote up, down and 
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strange quarks, respectively. The first term consists of the amphtude derived in the 
factorization approximation with some QCD corrections and the second one is a 
long-distance amplitude with c-quark or with M-quark in loop. At the hadronic level 
the second amplitude with the c-quark in loop can be interpreted in terms of the 
intermediate oi*^ D^*^ states which are frequently produced in _B-meson decays (see 
Fig. 1). 




Fig. 1. Example of Di*^D^*^ contribution to B decays into ttitK and KKK 



The charming penguin contribution to the B ■k'kK decay amplitude with two 
pions in P-wave is parametrized by the following expression: 

penguin — 2 G pTflpCpV p-j^-j^i^JTl-^.^ )pTrPKCOS0, (1) 

where Gp is the Fermi constant, rup is the p(770) mass, rp^Trl^inr) is the pirn vertex 
function of the tttt effective mass m,r7r , Pir and pk are the pion and kaon momenta 
in the p rest frame and 6 is the angle between the direction of flight of the tt^ and 
the direction of the tt+tt^ system in the B rest frame. The constant amplitude Cp 
reads: 

Cp = Ik A^-P{ml) {VubV:,Pu + VtbV*Pt) . (2) 

Here V^s are the Cabibbo-Kobayashi-Maskawa matrix elements, fx is the kaon 
decay constant, Aq ^''{m\) is the B ^ p transition form factor and Pu and Pt are 
penguin P-wave complex parameters to be fitted to the experimental data. 

Final state interactions in the isospin zero S'-wave are treated using the unitary 
model of Ref. ^^in which the tttt and KK channels are coupled. In this approach the 
sum of the several Breit-Wigner terms, usually used in the experimental analyses of 
the Dalitz plot distributions, is replaced by a set of unitary coupled meson-meson 
amplitudes. These amplitudes are expressed in terms of phase shifts Stt-k, Skk and 
inelasticity rj known from other experiments The P-wave pion-pion amplitude is 
well described by a Breit-Wigner term. No arbitrary phases nor relative intensity 
free parameters for the different resonances are needed. All the resonances appear 
in a natural way as poles of the meson-meson amplitudes. The scalar resonances 
/o(600) and /o(980) are examples of such poles of a single amplitude. The presence 
of the /o(980) is a dominant feature of the experimental m^rTr distribution. 
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3. Comparison with the experimental data 

We perform a fit to the data of Belle and BaBar^l Collaborations. Both groups 
have measured differential tttt effective mass distributions, branching fractions, di- 
rect CP violating asymmetries in charged B decays and time-dependent CP vio- 
lating asymmetry parameters in B'^ decays. Moreover, the Belle Collaboration has 
given helicity angle distributions in the p(770) and /o(980) regions. We obtain a 
good agreement with the data using only the four complex penguin parameters: 
Su = (0.15 ± 0.10) exp[i(1.90 ± 0.71)], St = (0.020 ± 0.002) exp[-i(0.26 ± 0.21)] 
for the S'-wave and P„ = (1.09 ± 0.21) exp[-i(0.98 ± 0.12)], Pt = (0.065 ± 
0.002) exp[-i(1.56 ± 0.08)] for the P-wave. 
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Fig. 2. The tt+tt effective-mass distributions in B^ — > tt+tt decays. The data points are 

taken from Ref. 4 and the sohd hnes represent the results of our model. 



In Fig. 2 we show a comparison of our model with the tt+tt^ effective mass dis- 
tributions measured by Belle for the B'^ ttttK^ and — > ttttK' decays. The 
large direct CP asymmetry of the order of 0.3, visible in the range of the p(770) 
resonance, is well described by our model. The fit is performed in the to^tt range 
between 0.60 GeV and 1.06 GeV where the two resonances p(770) and /o(980) domi- 
nate the pion-pion mass spectrum. Existence of these resonances leads to interesting 
interference phenomena seen in Fig. 3 where helicity-angle distributions are plotted 
for the combined B^ — > ttttK^ decays. In the p range the general behaviour of the 
data follows the cos^ 6 function, characteristic of a polarized p decay into tt+tt^. 
However, in the /o(980) range the angular dependence is not flat as one can expect 
for the decay of a S-wave resonance. It has an interference component proportional 
to cos 6 which originates from the presence of the p resonance tail under the /o(980) 
peak. 
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Fig. 3. Angular distributions in — » tt+tt" -K± (a) in the p(770) mass region and (b) in the 
/o(980) region. The data points are from Ref. 4 and the solid lines denote our model. 



We have also performed calculations of different observables for the B decays 
into three kaons in which a KK pair is in relative S'-wave. No extra free parameters 
have been used since all the parameters have been fixed in the fit to the B ttttK 
decays as described above. The averaged branching ratio for the B^ — > K'^K~K^ 
S-w&ve channels, integrated over the KK effective masses from the threshold up to 
1.1 GeV, is equal to 1.7x 10~^ and the corresponding value for the B° K+K~K° 
is 0.9 X 10~^. The direct CP asymmetry for the charged B decays is 0.07. The 
time-dependent asymmetry parameters for the neutral B decays are: S = —0.80 
and A = -0.13. 
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